Two-dimensional spin-1 frustrated Heisenberg model with valence-bond ground states 
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In this paper, we propose a method to understand the nature for the quantum disorder phase of 
the two-dimensional (2D) high spin frustrated modeL The ground state and excitation properties 
of a fuUy frustrated 2D spin-1 model are studied based on a model whose groundstate can be found 
exactly. By analogy to the pseudo-potential approach in the fractional quantum Hall effect, we 
conclude that the ground states of the fully frustrated spin-1 model are doubly degenerate valance 
bond solid (VBS) states along the horizontal or vertical direction of the square lattice. We also find 
that the ground state could be characterized by a nonzero string order, which rarely happened in the 
2D case. The method that we used reveals the connection between the fractional quantum hall effect 
and the frustrated 2D antiferromagnetism system. The VBS states capture the main character of 
the disordered phase in the 2D spin-1 frustrated system, and can be verified by a numerical method. 

PACS numbers: 75.10.Jm, 71.27.+a, 75.10.-b 



I. INTRODUCTION 

The 2D quantum frustrated magnets have attracted 
considerable attention in the past yearsi because they are 
believed to be promising candidates for realizing exotic 
spin liquid states. One of the most studied 2D frustrated 
models is the antiferromagnetic Heisenberg model on a 
square lattice with the first and second nearest neighbor 
interactions Ji and J2, 
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where (i, j) and ((i, j)) represent pairs of the nearest and 
next-nearest neighbors respectively. The general ground 
state properties of this model have been investigated in 
various methods^i^i^i^!^. It has been well accepted that 
the Ji — J2 model exhibits a quantum phase transition 
from a magnetically ordered Neel phase at small J2/J1 
to a spin disordered phase in the intermediate parameter 
region. Despite the success in predicting the order to 
disorder transition, the theory based on mapping onto 
the non-linear a model and symmetry analysis provides 
no clear information about the nature of the disordered 
phase. For the 2D spin-1/2 J1—J2 model, there have been 
a lot of works focusing on its nonmagnetic phase, but no 
definite conclusion has been drawn on the controversial 
question whether the disordered phase is a spin liquid 
without symmetry breaking or a valence bond crystal 
with broken spatial symmetry^'^. 

Most researches in this field focus on the spin-1/2 sys- 
tem because of its possible correlation with the High- 
Tc superconductivity. However, in recent years ultracold 
atomic system has provided an ideal playground to ex- 
perimentally investigate the high-spin strongly correlated 
system. Several proposals have been provided to realize 
the spin-1 lattice model-. Considering the rapid devel- 
opment in this field, we expect that it is possible for the 
spin-1 Ji — J2 model to be realized experimentally in 



the future. However, in comparison with the spin-1/2 
system, little knowledge is known about the disordered 
phase of the 2D spin-1 system except the possible dou- 
ble degeneracy of the ground state predicted by the field 
theorj*^. 

In this paper we first review some general results of 
the disordered phase in the frustrated model, then we in- 
vestigate the ground state properties of a fully frustrated 
spin-1 Ji — J2 model with J2/J1 = 0.5 by the pseudo- 
potential approach. Our method can be generalized to 
deal with other systems such as the frustrated honeycomb 
model and system with higher spin. 

The general property for the 2D spin 1 frustrated 
model was investigated in the scope of the (2-fl)- 
dimensional nonlinear a model"^'^, which is considered 
as the continuum limit of the SU(N) antiferromagneli^^. 
Though this model was derived in the semiclassical (large 
S) limit and no frustrated condition, it is still valid to de- 
scribe the condition of small S and frustrated case^.The 
effective action of a 2D spin-1 Ji — J2 model is 



cff 



1 



drd"r<j i(a^m)2 + c(Vm)^ }+Sb (2) 



and 



where Sb = "i-S Y.rVR § d'^R. ■ Ar, go 

c = ^8 (1 — 2a)JiSa with a — J2/J1. This is precisely 
the 0(3) non-linear cr model in the (2-1-1) dimensions 
with a residual Berry phase (Sb)- Via this model, renor- 
malization group (RG) analysis predicted that there is a 
transition from the Neel state to a disordered state and 
the transition point for the spin-1 model is ac = 0.46, 
larger than that in spin-1/2 case. The Berry phase plays 
a crucial role in determining the degeneracy and sym- 
metry of the ground state in the disordered phase^i^. 
For the spin-1 model, the degeneracy in the disordered 
phase is 2 and there is a spontaneous symmetry break- 
ing. We will see below that the symmetry of our trial 
ground states completely agrees with that predicted by 
the general theory above. 
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The analogy between the Heisenberg model and the 
fractional quantum Hall effect (FQHE) was first intro- 
duced by Arovas et al,^^ in the one-dimensional (fD) 
case. They decomposed the spin-f antifcrromagnetic 
Heisenberg Hamiltonian as: 

Hi=Y, Si-Si+i = ^[3P2(*, i+l)+P\t, i+l)-2] (3) 

i i 

while they treat the projection operator i + I) over 
spin-1 states associated with two consecutive sites as 
a perturbation and consider the Affleck-Kennedy-Lieb- 
Tasaki (AKLT) statei^, which is the exact ground state 
of the projection operator 3P^{i,i + l) over spin-2 states, 
as a trial ground state of the spin-1 Heisenberg chain. 
The AKLT state and the real ground state of the spin-1 
Heisenberg chain share a lot of important properties, for 
example both of them have an energy gap and the two 
point correlation functions decay exponentially. They are 
in the same universal class and characterized by a hidden 
order parameter called "string parameter" Numerical 
results show that the difference between the ground state 
energies of these two states is within 5%. In this work, we 
extended their method to the 2D fully frustrated model 
to study the properties of the groundstate and the ele- 
mentary excitation. 

For the fully frustrated model with J2/J1 — 0.5, we 
can rewrite the Hamiltonian Eq.(l) as: 

H = J2Ha with He, ^ J2 J2 ' (4) 

where the index a denotes the site of the dual lattice cor- 
responding to the ath plaquette. Each bond interaction 
between the two spins in a plaquette is equal. For con- 
venience, we expand this Hamiltonian by the projection 
operators of total spin in a plaquette: 

4 

H^/J2 = ^sPl ~ 4 (5) 
s=o 

with Cs = S{S + l)/2 and further decompose it as the 
summation of two parts: 

Ha = + H^, 

with 

H"JJ2 = 10P^-h6P^-4, 

The operator projects the spin state of the ath pla- 
quette onto the subspace with total spin S. One of the 
reasons for dividing the original Hamiltonian into two 
parts is that the model = J2a quasi-exactly 
solvable in the sense that its ground state can be obtained 
exactly. Furthermore, we observe that the coefficient Cs 
decreases rapidly as S descends. Therefore it is reason- 
able to apply the pseudo-potential method originally used 



in the FQHEis. Using this method, we first classify all 
the ground states of and treat = as a per- 

turbation. The ground state of is shown to be highly 
degenerate and belongs to two kinds of states: the dimer 
states and the valence bond solid states. However the 
huge degeneracy of the ground states is lifted when 
is introduced. Finally our result shows that the ground 
states of the fully frustrated spin-1 model are doubly de- 
generate and approximately described by the decoupled 
VBS states along the horizontal or vertical direction of 
the square lattice. We believe that our work would in- 
spire the interest in the spin-1 frustrated model. 



II. THE SPIN-1 MODEL WITH EXACT 
GROUND STATES 

Now we start by considering a spin-1 Hamiltonian on 
a square lattice: 

Hp = ^(c4P^+C3P^), (6) 
a 

where Hp is an arbitrary linear composition of projection 
operators and P^ with 04,03 > 0. Obviously, is 
just a special case of the above Hamiltonian. Since Hp is 
positive semidefinite, the state with the total spin of each 
plaquette < 2 is the exact ground states for Hp. Such 
a condition could be satisfied by two kinds of states: the 
dimer states and the valence bond solid states. 
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FIG. 1: (Color online) Different dimer ground states of Hp. 
The green (dark) ellipsoids represent singlet dimers formed 
by two S = 1 spins. 



The dimer states are the states with each plaquette 
possessing a dimer, which is a spin singlet formed by two 
neighboring spins: ^(|1, -1) + | - 1, 1) - |00)). This 
family of states are illustrated in Fig.l and are very sim- 
ilar to the dimer gromid states in the spin- 1/2 case^ 
except the singlet dimers are formed by pairs of spin 
with S = 1 and also the defect spin in Fig. 1(d) has 
S = 1. As stated by Batista et ai, the states in Fig.l (a- 
d) represent all the possible configurations which satisfy 
the constrain that every plaquette possesses at least one 
dimer—. For the configuration in Fig. 1(a), rotating any 
array of dimers along a diagonal direction by 7r/2 , we 
will get a new degenerate state of this model, as stated 
ii>i^. The degeneracy of this kind of ground state is pro- 
portional to 2^^+^, where N is the total number of the 
sites. States in Fig.l(b-d) represent the configurations 
with one "defect" and it is possible to have a localized 
S = 1 spin as illustrated in Fig. 1(d). Translating or 
rotating these kinds of states around the central defect 
gives a huge amount of degeneracy proportional to N, the 
total number of the sites. Configurations in Fig. l(a-d) 
and their corresponding degeneracy states provide a full 
configuration subspaces that satisfy the constraint that 
each plaquette has at least one singlet. 

For the spin-1 system, however, there may be another 
family of ground state totally different from that of the 
spin-1/2 case: the VBS state, or the AKLT state^^. An 
S = 1 spin in one site could be considered as a symmetric 
composition of two = 1/2 spins and the AKLT state 
is formed by the valence bond, which is the singlet state 
of two neighboring S" = 1/2 spins. If there is a global 
state ensuring that each plaquette has at least two va- 
lence bonds, as illustrated in Fig.2(a-f), such a state is 
also a ground state of Hp due to S*^ < 2. Obviously, this 
kind of ground state is also highly degenerate. The lat- 
tice we considered is under the open boundary condition 
with the total site of lattice N oo. Configurations in 
Fig.2(a-c) are composed by many ID AKLT chains with 
infinite length. Fig.2(d-e) represent the configurations 
composed by the closed AKLT loop and a defect in the 
center. Other configurations not listed here which satisfy 
the same constraint are similar to Fig.2(e,f), but with a 
defect of a different shape. We could estimate that the 
number of this kind of degeneracy is proportional to N 
by considering the possible positions of the central defect 
in Fig.2(d-f). 

Now we focus on the excitation of the second fam- 
ily. The elementary excitation of the AKLT chain is the 
magnon excitation with a Haldane gap of 0.35^"'^, in which 
a singlet dimer formed by two neighboring S* = 1 /2 spins 
is excited to a triplet state, as illustrated in Fig. 3(a). 
Unlike the dimer case, these two 1/2 spinous could not 
be separated without increasing the energy proportional 
to the distance between them, so this kind of excitation 
is confined. However, the triplet state itself, which is 
formed by two binding spinous, could propagate along 
the AKLT chain freely, as shown in Fig. 3(b). 

It is noted that a similar method has been used to con- 
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FIG. 2: (Color online) Different AKLT ground states of Hp. 
The blue (dark) ellipsoids represent singlet states formed by 
two S — 1/2 pseudospins. 



struct the exact ground state of the well known ID AKLT 
modeli^ and the 2D spin-1/2 model with additional cyclic 
exchangesi^. The application of such a method has made 
a lot of success in understanding the nature for spin 
systemsiS^^^il^. 



III. THE GROUND STATE OF THE A - Ja 
MODEL WITH FULL FRUSTRATION 

Since the Hamiltonian Eq.([6]) is somehow an artificial 
model, how could we expect that it is related to the fully 
frustrated Ji — J2 model? We will clarify this problem 
below using the analogy to the pseudopotential approach 
in the FQHE^'t, where, restricted within the first Lan- 
dau level, the interaction can be expanded by relative 
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FIG. 3: (Color online) (a) One singlet is excited to a triplet 
(open ellipsoid), (b) The excitation could propagate along 
the AKLT chain freely. The blue (dark) ellipsoids represent 
singlets formed by two pseudospins with S — 1/2. 



angular momentum projection operators: f (r^ — rj) — 
J2'k'=Q'^kPkiij) with Pkiij) the projection operator on 
states with angular momentum k. If a potential satisfies 
Vk=0 when k > m, the LaughUn state is the exact ground 
state for this model. The success of this method depends 
on the difference between the coefficient Vk of fc = m — 2 
and k = mi^. 

Now we treat the 2D Ji — J2 Hamiltonian in the spirit 
of the pseudopotential approach and consider as per- 
turbation of just as Arovas et. al did in the ID case. 
Since the ground state of = J2a been ana- 

lyzed above, all the degenerate ground states of Hamil- 
tonian ^ may be a possible candidate of the trial ground 
state for the Hamiltonian However, the perturbation 
= i?^ lifts the high degeneracy and chooses some 
states as the best trial ground state. 

Next we will calculate the energy expectation under 
those degenerate states. The lower the expectation value 
is, the closer the state is to the real ground state of 
We notice that our lattice is constructed by the bonds 
connecting one site with its nearest neighborhood or the 
next nearest ones. To calculate the average energy of 
various candidate trial groundstates, we must first calcu- 
late the energy expectation of each bond. For the dimer 
states as shown in Fig.l, all the configurations in Fig.l 
(a-d) are constructed by three kinds of bonds. Take the 
Fig.l (d) for an example: the first kind is represented by 
bond 12, which is covered by a single dimer. Bond 23 
represents the bond between two different dimers, and 
the third kind (bond24) is the one between a dimer state 
and a single spin with S ~ 1. The energy expectations 
for these three kinds of bonds are 



(^1 ■ ^^)dv 
(^2 ■ ^s)dv 



(7) 



respectively. It is ready to calculate the average energy 
per plaquette in the thermodynamic limit: 



{E) 



dimer 



= -2J2, 



which is the same for all the configurations in the Fig.l 
(a-d). 



Now we turn to the VBS states in Fig. 2. Similar to 
the method above, we can classify all the bonds which 
construct the configurations in Fig. 2 into five classes: (1) 
Bond 12 or 34 in Fig. 2(a), which is covered by a va- 
lence bond formed by the singlet of two 1/2 spinous in a 
AKLT chain with infinite length; (2) Bond 57 in Fig. 2(c) 
is the bond connecting two next nearest sites in an AKLT 
chain; (3) Bond 13 or 24 in Fig. 2(a) which connects dif- 
ferent AKLT chains; (4) Bond 89 in Fig. 2(d) is similar to 
the first one, but the bond is covered by a valence bond 
in a closed AKLT loop with finite length; (5) Bond 90 in 
Fig. 2(d) is the bond between one site in an AKLT chain 
and another site occupied by a single spin with S = 1. 

Now we calculate the average energy under all these 
candidate trial ground states in Fig. 2 to find which one 
has the lowest energy. We define as an operator 

to rotate the system by 7r/2 around one site in the lat- 



tice. So we can see that ^^4^^2-11 



--T^,2'^\klt\\ ' and 



the Hamiltonian in Eq.(l) is invariant under the operator 
So the states in Fig. 2 (a) and (b) are degenerate. 
We denote these states as aklt\\- In a similar way, we 
can get 



(S, -S 



^lAKLTW 



^ilAKLTW 



-4/3 



while other bonds all contribute zero to the average en- 
ergy. It is clear that we have 



{E) 



AKLTW 



which is smaller than (E) dimer- If we consider the states 
like Fig. 2(c), where there are many right angles along 
the AKLT chain, we find that the average energy under 
these states denoted by '^aklt± is always higher than 
that of AKLT\\- In the AKLT state, there is a short 
range correlationi^, and thus we get 

(S, . S,) = 4/9, 

whereas we have (Sg -Sg) = (Sg-S^) — —4/3. So this kind 
of bonds in the corner of the AKLT chain would always 
raise the average energy of the system. We can find in 
Fig. 2(d) (Sg • Sq) = So the energy of configurations 
corresponding to Fig.2(d-f) and other configurations not 
listed here is always higher than the energy of ^^^iry 

and '^''aju^j'w for two reasons: the appearance of right 
angles along an AKLT chain as well as the fact that the 
average value of Sj • Sj+i is lower in an AKLT chain with 
infinite length than that in a closed AKLT loop with 
finite length, i.e., {S^-S^)aklt\\ < {^s'^^) akltu- The 
proof of this fact has been shown in the Appendix. 

As we showed above, the perturbation lifts the degen- 
eracy and chooses state ^^^^j^n and ^^^^^^n (shown 
in Fig. 2 (a, b)) as trial ground states. Apparently, they 
are doubly degenerate. It is natural to ask how close 
they are to the real ground states. Here we give an es- 
timation. We can expand ^^^^^t^h in the spirit of the 



pseudo-potential approach: 



AKLTW 



^a^2 + ^^1, where 
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= Ec^f^Wll (/5 = 1-2) and + 62 = i. The 
reason why there is no component of is that if there 
exists a state in which the total spin of every plaquette 
is zero, the energy for every plaquette gets to its min- 
imum simultaneously which is impossible because the 
wave function would not be self-consistent. Observing 
that {^'iotai) AKLT\\ = &a -\- 2h = 8/3, we can get the co- 
efficient b = 0.994, which means (*1klt||I*i> = 0-994. 
Based on this result, we can reconfirm the validity of the 
perturbation used above. Because ^'i is the dominant 
part in '^\jii^rp^\^, the operator in the perturbation 

plays a more important role than P^. Furthermore 
since the coefficient of P^ is smaller than that of P^, 
the effect of the perturbation to the unperturbed state is 
even smaller. Therefore, we believe that this trial state 
should be a good approximation to the ground state of 
Hamiltonian 

Nevertheless, it seems that there is still a question to 



be clarified. Since the trial ground states VE"^ 



AKLT\ 



and 



^ AKLT\\ ^'^^ degenerate, we should use the degenerate 
perturbation rather than that used above, which means 
that a linear composition of these two degenerate states 
could further lower the energy and forms a new ground 
state, just as in the resonant valence bond theory. We 
shall show that it is not true in this case. As we know, 
whether the degenerate perturbation works depends on 
the non-diagonal term {^\xlt\\\-^\'^\klt\\) vanishing 
or not. For an iV x square lattice model with periodical 
boundary condition, the wavefunction corresponding to 
Fig. 2(a) can be represented as 



^ AKLTW 



^ AKLTW 



^ AKLTW 



^ AKLTW ' 



(8) 



with ^"aklt\\ denoting the wavefunction of the 
71— th periodical fD AKLT chain. We estimate 
^'^°aklt\\\^\'^\klt\\) ^ N'^2~^\ which decreases ex- 
tremely rapidly when N increases. So we can safely draw 
the conclusion that in the thermodynamic limit, the non- 
diagonal term is zero and the degenerate perturbation 
and non-degenerate perturbation give the same result. 



IV. EXCITATION 

We can also get some results of the low energy excita- 
tion based on our trial ground state wave function. Taken 
the '^'Xkltw ^'^^ example, the basic picture of the low en- 
ergy excitation is illustrated by Fig. 3 where a singlet in 
an AKLT chain is excited to a triplet and propagates 
along the chain. Now we calculate some quantitative re- 
sults of this picture using a variational technique. Since 
our trial ground state wave function could be written as 
the production of many parallel ID AKLT chains and 
the correlation function between different AKLT chains 
is zero, the low energy excitation is similar to the ID 
caseiiii^. To find the variational wave function of the low 
energy excitation, we use the matrix product statCf^S, At 



a site n, we define the matrix 



/T73|0)„ y^\-l)„ 

/2/3\l)n -7173 |0)„ 



(9) 



In terms of the matrix product, the ground state of the 
i—th AKLT chain could be represented as 



N 



1^ AKLT\ 



(10) 



and the overall ground state is given by Eq.®. 

The excited state as shown in Fig. 3(a) could be repre- 
sented as 



' AKLT\\ ■••'^ylKLTII ("-)•■■ ® '^AKLTW 

(11) 

where ^^'^^^t^h (n) denotes the state with the n— th sin- 
glet of the i—th AKLT chain being excited to a triplet. 
For brevity, we write \^'^xlt\\(''^^'^ I"") '^hich can be 
represented in terms of the matrix product as 



l/3|0)„+i v/273|-l)„+i 
2/3 |l)„+i yr/S |0)„+i 




(12) 



m—n+2 



This is a state with Sf^^. = 0. After some algebra, we find 
that 

{m\n) = (5m,n + -((5„i,n-l + (^m.ri+l) (13) 



(mlHln) 



8Ji^ 
~9~' 



Ji 



) (14) 



where H is the Hamiltonian (4) or (1) with Ji — 2J2. A 
momentum eigenstate is thus defined as 



n 

which satisfies 



-AKLTl 



(n)), 



(k|k) = {l + -cosk,)N\ 



2 

— ( 

'9 9 



(k\H\k) = {- - -cosk^)N^Ji. 



So the variational energy is 



8 - 2cosk^ 
^-^^^ - 9 + 2cosk/'^ 



(15) 

(16) 
(17) 

(18) 



where we have made an energy shift of —^J2N'^ in the 
above calculation. 
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V. DISCUSSION AND GENERALIZATION 

The string order was first observed in the spin-l AKLT 
chain^^. It is behaved that even when the Hamihonian 
deviates from the exact AKLT point, the string order did 
not vanish^i. However, the nonzero string correlator was 
rarely observed in the 2D model, the only exception that 
we have known up to now is the Wen-Kitaev Mode^, 
which is a 2D exactly solvable model with a string cor- 
relator attaining finite value. In our trail ground state, 
since there is a strong evidence that the 2D ground state 
is decoupled into several ID AKLT chains, the order 
parameter— proposed to characterize the ID VBS state 
should also be valid in our 2D problem. We believe that 
the string order parameter: 

(19) 

should get its maximum in the point J2/ Ji — \I2, which 
is still to be verified by numerical methods. 

Now we discuss the case when the coupling coefficient 
J2/J1 deviates from 1/2. For the general Ji — J2 model, 
the RG analysis based on the nonlinear a model can show 
that J2/J1 — 1/2 represents a stable fixed point in RG 
flow in the parameter spaced, so we believe that our trial 
ground states represent a universal class and capture the 
basic property of the spin disordered phase of the spin-l 
Ji — J2 model. 

VI. CONCLUSION 

Based on a solvable spin-l model with dimer-type and 
valence-bond-solid-type ground states, we studied the 
fully frustrated spin-l Ji — J2 model and proposed a pos- 
sible trial ground state at the point J2/J1 = 0.5, which 
has been seldom studied before. Such ground states are 
doubly degenerate and approximately described by the 
decoupled AKLT states along the horizontal or vertical 
direction of the square lattice. We believe that this trial 
state captures the main character of the disordered phase 
in the 2D spin-l frustrated system, and can be detected 
by numerical methods. In addition, the pseudo-potential 
method used here is not only restricted to square lattice 
and the spin-l case, but also to higher spin or other lat- 
tice cases as long as the Hamiltonian could be expanded 
to the summation of projection operators. Taking the 
spin-l Honeycomb lattice for an instance, the state with 
gtotai ^ 2 for each hexagon is close to the real ground 
state of the Ji — J2 — Jd, model in the Honeycomb lattice, 
even closer than the case of the square lattice because the 
perturbation takes a smaller portion in this Hamiltonian 



than that in the square lattice Ji — J2 model. So our 
method offers a potential framework to explore quantum 
exotic state in spin systems. 

This work is supported in part by NSF of China under 
Grant No. 10574150 and programs of Chinese Academy 
of Sciences. 

APPENDIX A 

In this appendix ,we will prove the fact that the average 
value of S. • S.-)_i is lower in a AKLT chain with infinite 
length than that in a closed AKLT loop with finite length. 
Our proof is based on some results in Reii^. 

Let (0),and (— ) denote the orthonormal basis for 
spin 1 consisting of eigenstates of S^' with eigenvalues 
-1-1,0, and —1, respectively and 

l/-!! = \/2( + ), ^i2=V21 = (0), ^22 = %/2(-). (Al) 

We denote the wavefunction of a single AKLT chain with 
length L as 

(A2) 

where e°'^ is the antisymmetric tensor with e^"^ = 1. As 
shown inJ^, the inner product of two AKLT chains with 
L sites is 

■ n'^ = <5:;<5f (3^ - 1)/2 + sp'^, (A3) 

so the normalization of the ground state with periodic 
boundary conditions reads 

nf • = 3^ + 3 (A4) 

Next we calculate the spin-spin correlation function 
between the site 1 and 2 (S'f S"!) for a chain with L sites 
and periodic boundary conditions. FoUowing^^, the spin 
operators acting on our basis for a single spin 1 give rise 
to 

S^Vaft = (l/2)[(-l)" + (-l)''^]^„^„ (A5) 
^2^a.ft = (-l/2)[(-l)"^ + (-l)^^]V^„,^,. (A6) 

Combining A{2) — A{6), we can find 

{S^Si) = -ii3^-'^-l)/{3^ + ^). (A7) 

We notice that the AKLT chain is isotropic, so {S^ ■ 
S^)=3{SfSi) = -12(3^^2 _ lyi^^L j^Q^ jjj^yg 

proven the fact that the average value of S( • S,+i is lower 
in an AKLT chain with infinite length than that in a 
closed AKLT loop with finite length. When L — > 00, we 
get the value in^^ ^i^a,t (S. • S.+i) = -4/3. 
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